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Abstract. We analyze a flux-limited sample of persistent and bright (with 2-10 keV fluxes exceeding 1.4 x 10~ 10 
erg s" 1 cm -2 ) low-mass X-ray binaries (LMXBs) in our Galaxy. It is demonstrated that the majority of binary 
systems with X-ray luminosities below log L(erg/sec) ~ 37.3 have unevolved secondary companions (except for 
those with white dwarf donors), while systems with higher X-ray luminosity predominantly harbor giant donors. 
Mass transfer in binary systems with giants significantly shortens their life time thus steepening the X-ray lumi- 
nosity function of LMXBs at high luminosity. We argue that this is the reason why the LMXB luminosity function 
constructed in the last years from observations of sources in our and distant galaxies demonstrates a break at 
logL(erg/sec) ~ 37.3. 



1. Introduction 

Binary systems with low-mass secondary companions to 
a compact star (black hole or neutron star) were discov- 
ered in the 1960s at the dawn of X-ray astronomy. All-sky 
surveys performed by different orbital X-ray observatories 
(UHURU, HEAOl, Ariel V, etc.) provided us with a rel- 
atively large sample of such objects in our Galaxy. The 
advent of focusing X-ray telescopes with an angular reso- 
lution of arcseconds initiated studies of such low-mass X- 
ray binaries (LMXBs), first in nearby galaxies, like M31 



(Trinchieri & Fabbiano 1991 Primini, Forman, & Jones 



19931, and then in more distant galaxies (see e.g. reviews 
Kim fc Fabbiano|[2004{|Gilfanov|[2004 |. 



The accretion luminosity L x of a persistent (sub- 
Eddington) LMXB is directly proportional to the mass 
transfer rate M from the secondary star. According to 
the standard theory of binary star evolution, the persis- 
tent mass transfer via the inner Lagrangian point in a close 
binary is due to an increase of the size of the donor star 
relative to its Roche lobe. In an LMXB, this can be done 
either by decreasing the size of the Roche lobe due to loss 



of the orbital angular momentum (see e.g. Paczynski & 
Sienkiewicz||1981| |Verbunt fc Zwaan||1981[ ) , or by increas- 
ing the radius of the star as a consequence of its nuclear 



evolution ( 


Webbink, Rappaport, & Savonije 


1983 


Taam 


1983 


Ritter 


1999 


). The observed properties of different 



populations of LMXBs can be used to test the models of 
binary evolution. 



The X-ray luminosity function (LF) is an important 
characteristic of the LMXB population. In other galax- 
ies it is directly constructed from measured X-ray fluxes 
down to a luminosity L x ~ 10 37 erg/s and can be fitted 
by a power law dN/dlogL oc £,—0-8— -1.2 w jth a steep- 
ening at luminosities exceeding the Eddington limit for 



accreting neutron stars (logL x > 38.5) (Kim & Fabbiano 
2004). However, recently it became possible to construct 



the luminosity functions of LMXBs in nearby galaxies, 
like M31 or Cen A, down t o much smaller luminosities 
of the order of 10 36 



erg/s (Primini, Forman, & Jones 
19931 |Gilfanov| [20041 |Voss fc Gilfanov||2006| |2007| |Voss 



et al. 2009). Over such a wide luminosity interval the lu- 
minosity function of LMXBs can no longer be described 
by a single power law and demonstrates a characteristic 
break at \ogL < 
down to luminosities L 



37.3 ( |Gilfanov||2004[ ). 

10 3tl erg/sec was obtained for 



A similar result 



galactic LMXBs from all-sky surveys (Grimm, Gilfanov 



& Sunyaev 2002 Gilfanov 20041 and the survey of the 



Galactic bulge (Revnivtsev et al. 2008). The statistical 



significance of the LF break in the above-mentioned works 
is high, therefore the existence of the break requires phys- 
ical explaination. One of the clearest case of the flattening 
of the LMXB LF at low luminosities can be seen in the 
work of Revnivtsev et al. (20081, where the LMXB can- 



didates in the bulge of our Galaxy were traced down to 
luminosities 10 35 erg/sec, almost unreachable for LMXBs 
in outer galaxies. 

Different explanations have been proposed for the ori- 
gin of the break in the observed LF of LMXBs. For 
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example, the break in the LF might be caused by the 
change of the dominant orbital angular momentum loss 
mechanism from magnetic stellar wind ( |Verbunt fe Zwaan 
1981 1 to gravitational wave emission (see e.g. |Paczynski 
& Sienkiewicz 1981 ) in the population of LMXBs with 



low-mass main sequence secondaries ( Postnov & Kuranov 



2005 1 . It was also noted already quite long ago that the 
most luminous LMXBs with high mass transfer rates 
(M > 10 -8 -Mo/year), could harbor giant donors (see 



e.g. |Webbink, Rappaport, fc Savonije|1983 1 , and this may 
undcrly the difference between the low- and high- luminos- 
ity sources. The study of the LMXB LF in galaxies using 
population synthesis methods could potentially be quite 



powerful (see e.g. Fragos et al. 2008 Kim et al. 2009) 



However, this approach involves a number of uncertain 
parameters of binary star evolution, so it is hard to make 
firm conclusions based on the population synthesis simu- 
lations only. 

In our paper we make an attempt to identify the main 
physical reason for the origin of the observed break of the 
broad band luminosity function of LMXBs. 

We show that the break in the LF of LMXBs is prob- 
ably caused by different types of donor stars: evolved 
secondaries (giants) at luminosities above the break at 
log L(erg/sec) ~ 37.3, and main sequence stars at lower 
luminosities. The correctness of our conclusion might be 
checked by direct calculations of the LMXB properties us- 
ing the methods of population synthesis, which we plan to 
do in the future. 

2. A Flux limited sample of persistent LMXBs in 
the Galaxy 

We consider only persistent LMXBs. This selection can 
be relatively easily done for galactic LMXBs, but in gen- 
eral this can be a difficult task for LMXBs in outer 
galaxies, since they are typically observed only during 
short time intervals. Monitoring of X-ray sources in our 
Galaxy shows that at any particular time there are only 



a few bright transients (see e.g. Remillard, Levine, & 
McClintock 20091, whose effect on the instantaneous LF 



of all galactic LMXBs is rather small. 

Moreover, the properties of the observed luminosity 
distribution of individual variable LMXBs were found to 
not affect the shape of the instantaneous LF of the galactic 
sample ( Postnov fc Kuranov||2005 ). Therefore the results 
of our present study can also be applied to the instanta- 
neous sample of LMXBs in distant galaxies. 

We study a representative sample of persistent galactic 
LMXBs with X-ray luminosities L x > 10 36 erg/s and com- 
pare their properties above and below the break. Our goal 
does not require the sample to be volume limited or com- 
plete, because we are not interested in statistics (which 
was extensively studied in works of other authors), but 
rather we look for a correspondence of the time averaged 
X-ray luminosity of sources to some parameters of their 
binaries. 



To select persistent Galactic sources we can analyze 
any existing X-ray survey. However, such a survey should 
fulfill the following requirements: 1) it should cover as 
much as possible of the Galactic plane region, where the 
majority of the LMXBs are located, 2) it should have an 
angular resolution of better than ~ 2 — 3° in order to avoid 
confusion within the Galactic plane, 3) it should be per- 
formed in the energy band ~ 1 — 20 keV where Galactic 
LMXBs emit most of their bolometric luminosity. These 
requirements significantly shorten the list of usable sur- 
veys. In particular, the sky survey of the RXTE obser- 
vatory ( Revnivtsev et aL]|2004 1 does not cover the region 
of the Galactic plane, and the scans of the RXTE ob- 
servatory over the Galactic bulge and the Galactic plane 



(Markwardt et al. 20001 have not (yet) been used to 
perform a sky survey, rather it is used to measure the 
fluxes of preselected list of sources. Surveys carried out 



by the ASCA (Sugizaki et al. 2001 Sakano et al. 20021 



and BeppoSAX (Sidoli, Belloni, fc Mereghetti 2001) ob 



servatories cover only a small fraction of the Galaxy. Even 
smaller sky areas are covered by the Chandra and XMM 
observatories. The survey of the INTEGRAL observatory 
( |Krivonos et al. [2007 ) covers the whole Galaxy, but it is 
done at hard X-ray energies, and is not fully suitable for 
our purposes (however, we will use it for additional checks 
of "persistency" of selected sources). 

So the persistent X-ray sources detected by UHURU 



(Forman et al. 1978) turn out to be most suitable for 



our study; as an additional check of their persistent 
behavior we have also examined their presence in the 



INTEGRAL all sky survey ( |Krivonos et al^[2007| ) . We se- 
lected only sources with 2-10 keV fluxes above 1.4 x 10 -10 
erg s _1 cm~ 2 which ensures that their luminosity is not 
lower than 2 — 3 x 10 36 erg/s up to distances of 12 kpc (i.e. 
up to the further edge of the Galactic bulge). This allows 
us to conclude that we do not miss brighter sources at 
least in more than a half of the Galaxy. We also included 
in our sample the black-hole binary GRS 1915+105, which 
was not seen by UHURU during its operation time (1971- 
1973), but now remains persistently bright since its ap- 
pearance in 1992 ( |Castro-Tirado, Brandt, fc Lundfl992"l ) , 
and excluded the source GX 1+4, which is known to 



be accreting via a stellar wind (Hinkle et al. 20061 and 
not via Roche lobe overflow which we consider here. 
The fluxes assigned to GRS 1915+105 and GS 1826-24 
were taken from measurements of the All Sky Monitor of 
the RXTE observatory averaged over 1996-2009 (see. e.g. 
http : //xte . mit . edu/ asmlc/ASM . html ) . 

Energy fluxes from the selected sources were calcu- 
lated from observed UHURU count rates, assuming that 
the Crab nebula count rate in the energy band 2-10 keV 
corresponds to the energy flux 2.22 x 10 -8 erg s _1 cm~ 2 . 

In Table 1 we present the list of sources with values of 
their orbital periods and estimates of their distances and 
luminosities. 
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Table 1. The brightest persistent low mass X-ray bina- 
ries in the Galaxy with fluxes > 1.4 x 10~ 10 erg s _1 cnr 2 



measured by UHURU (Forman et al. 1978). Orbital pe- 



riods of systems are adopted from Ritter & Kolb (2003) 



and Liu, van Paradijs, & van den Heuvel (20071 



Name 


L x (2-10 keV) 
10 37 erg/sec 


Dist, kpc 


Period, h 


Sco X-l 


37.4 


2.8 1 


18.94 


GX 5-1 


26.1 


9 2 




GRS 1915+105 


26.0 


ll 3 


739.20 


Cyg X-2 


19.4 


11. 6 5 


236.27 


GX340+0 


18.2 


12 4 




GX 349+2 


15.3 


8.5 6 


22.5 


GX 17+2 


14.9 


7.5 7 




GX 9+1 


8.7 


7.2 7 




4U1820-30 


5.7 


8 8 


0.19 


Ser X-l 


5.5 


8.4 9 


12.96 


GX 13+1 


5.5 


6.9 10 


601.7 


4U 1735-44 


4.9 


9.1 11 


4.65 


GX 3+1 


3.4 


4.5 12 




4U1624-49 


3.1 


15 13 


20.9 


4U1636-53 


2.4 


5.9 10 


3.79 


GX 9+9 


1.8 


5.0 7 


4.20 


4U 1746-37 


1.4 


11.0 s 


5.16 


4U 1705-32 


1.18 


13 14 




1A1742-294 


0.81 


8.5 15 




4U1254-69 


0.79 


ll 16 


3.93 


4U0513-40 


0.74 


12.1 s 


0.28 


4U1823-00 


0.61 


6.3 17 


3.19 


4U1915-05 


0.43 


8.8 s 


0.83 


4U0614+09 


0.34 


3.2 18 


0.81 


4U1702-42 


0.32 


6.2 19 




4U1626-67 


0.32 


g 20 


0.69 


GS1826-24 


0.31 


6.0 21 


2.25 22 


4U 1708-40 


0.27 


g 15 




4U1543-62 


0.26 


7.0 s 


0.30 


4U1724-30 


0.17 


9.5 s 




4U1850-08 


0.17 


8.2 s 


0.34 


4U1812-12 


0.09 


4.0 19 




4U1556-60 


0.09 


4.0 7 




4U1822-37 


0.04* 


2.5 23 


5.57 
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Fig. 1. Positions of the LMXBs from Table 1 on the or- 
bital period - X-ray luminosity diagram. Filled circles de- 
note the positions of peculiar systems: gray circles show 
the positions of binaries with white dwarf or hydrogen- 
deficit accretors, the dashed circle shows the position of 
4U1822-37 in which we see only a fraction of the total 
X-ray luminosity due to its nearly edge-on orientation (a 
so-called accretion disk corona source). The region, oc- 
cupied by binary systems with a giant donor is shown 
by the shaded squar e. The dashed line sho ws the slope 
L x oc P 2 ' 5 , derived byllben & Tutukov|(|l984| for the mass 



tranfer rate in binaries in which angular momentum loss 
is driven by a magnetic stellar wind from a main sequence 
donor. The dotted line shows the boundary of persistency 
according to the thermal-viscous disk instability for bina- 
ries with low-mass main sequence donors L x /1 37 erg/sec 
- 0.025 P hr 4 from eqn.32 of |Dubus et al] ( |l999[ ), where we 
adopted Mi = 1.4M , M 2 = O.4M and L x = O.lMc 2 . 



(1) -|Bradshaw, Fomalont, fc GeldzahTer||1999|(2) - 1 Jonker et al. 
1 2000 1 (3) -IFender et al.||1999| (4) - I van Paradijs fc White||1995| 
(5) -|Smal e |l55§ | {ti) -|Wachter fe MargonllWI (7) -Christian & 
Swank|1957[ ( 8) - fKuulkers et al.|2U03| (9) -Irjbisirzaki, Sug imoto, 
Hllm^rwa^^g^ (10) — ^ancTyopaTin^ ™ 
et al.Hitol (12) - 



Cuulkcrs &; van der 



I.fi9^| (1 1) - |Augustel.jn' 
kUs||2000| (13) -|xfang7 



fc Nowak|2007| (14) - |in't Zand, Cornelisse, Mendez|i|2'oogt 
(15) - adopt ing dist ance to the Galactic Center, (lb) — |(Jourvoisier| 
|et al.| <1986), (17) -|Shahbaz, Watson, fe Hernandez- Feralta||2007j 
(lij) Ifkuulkers et af. | liOi(J| (It)) - 1 Jonker k Nelemans|!i0u4| (20) - 
|(Jriakrabarty|1998| (21) - i Heger et al.|20U7| (22) - Mescheryakov et 
al. 2010, (23)- |Mason fc Cordova|1982| 

* — the source belongs to a group ot so called accretion disk corona 
sources, in which wc do not sec the direct emission of the innermost 
parts of the accretion flow, therefore the observed X-ray luminosity 
of such a source is only a part of its intrinsic X-ray luminosity. 



3. Results 



apparently have large sizes. This (in combination with 
the assumption that the companion star fills its Roche 
lobe) means that donor stars in these systems must have 
large radii and thus are evolved stars (subgiants or gi- 
ants). In several cases this was shown directly from the 
observed spectra of the optical counterparts (Sco X-l, 
|Bandyopadhyay et al.||1997| GX 349+2, [Bandyopadhyay 
|et al. | |1999[ Cyg X-2, |Casares, Charles, fc Kuulkers||1998 ' 



GX 5-l, |Bandyopadhyay, Shahbaz, fc Charles||2003| GRS 
1915+105, |Greiner et al.|2001[ Cir X-l, | Jonker, Nelemans,' 
fc Bassa||2007[ ). 

To the contrary, there are virtually no persistent 
LMXBs with X-ray luminosities lower than ~ 10 37 5 erg/s 
and orbital periods longer than 5-10 hours (see also|Ritter 



From Table 1 and Fig.l it is clearly seen that all systems fc Kolb[2003 Liu, van Paradijs, & van den Heuvel 
with X-ray luminosities exceeding ~ 5 x 10 37 erg/s for 
which we have information about their orbital periods, 



2007). 



This is consistent with the assumption that donor stars 
in these systems are main sequence stars (or degenerate 
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dwarfs, if the orbital period is smaller than ~ 1 hour). 
For cataclysmic variables with accreting white dwarfs and 
orbital periods less than ~10 hours (which are much 
more numerous in the Galaxy and can be found closer 
to the Sun, so we can study them in much greater details) 
this statement has been solidly supported by observations 
( Smith fc Dhillon||1998[ ). Therefore, we conclude fi •om ob- 
servations that at luminosities L x ~ 10 375 erg/s there is a 
transition from LMXBs with predominantly giant donors 
to binaries with predominantly main-sequence donors. 

4. A simple model of the LMXB luminosity 
function 

Numerous studies have been made in the past to describe 
the evolution of the LMXB population using specific mod- 
els of LMXB evolution. They were mostly aimed at un- 
derstanding the evolution of the total number of galactic 
LMXBs without distinguishing their X-ray luminosities 
(see e.g |White fc Ghosh] p98) [Ghosh & Whitel |200ll> 



Attempts to describe the LMXB LF can be found in 



Webbink, Rappaport, & Savonije (19831, but this work 



concentrated only on the bright end of the LF and does 
not consider the effects of magnetic braking on the pe- 
riod/luminosity distribution of binaries. The importance 
of the magnetic braking for relatively short period sys- 
tems was emphasized, for example, in papers by |Pylyser| 
fc Savonije|(|1988|);|Cote fc Pylyser|(|1989[). In these papers, 



an indication of transition from a population of bright sys- 
tems dominated by nuclear evolution of donor stars to a 
population of systems in which mass transfer is sustained 
by the magnetic braking can be found. At that time it 
was impossible to make a quantitative comparison of the 
model predictions with the observed samples of LMXBs 
due to the lack of sensitive galactic X-ray surveys. So 
these studies were mostly focused on the exploration of 
evolutionary tracks leading to the formation of particu- 
lar objects, like rotation-powered pulsars (e.g. |Kulkarni 



fc Narayan|p88| |Naylor fc Podsiadlowski||1993| |Deloye 



2008 and references therein), or ultracompact binaries 



2002). 



(e.g Nelemans et al. 2001 Belczynski, Kalogera, & Bulik 



Many papers that appeared later on were de- 
voted to the population synthesis of LMXBs (see e.g. 



Belczynski, Kalogera, fc Bulik|2002 Pfahl, Rappaport ,~fc" 



Podsiadlowski||2003[ |van der Sluys, Verbunt, fc Pols|2005 



Fragos et al.|2008 ), which until now suffer from uncertain- 
ties in the binary star evolution (e.g. the treatment of the 
common envelope stage, parameters of magnetic braking, 
etc.). Clearly, such detailed studies are very important and 
must be continued as they may help to constrain values of 
specific parameters of binary evolution from comparison 
with observations. However, to understand the observed 
gross properties of the LMXB LF we can try to use very 
general assumptions about binary star evolution, as we 
show below. 

Consider an ensemble of binaries that has been formed 
in a galactic halo over a time period T in the past. Let 



t(P) be the duration of an active LMXB stage of a source 
within the orbital period range P, P+dP. Let n(P) be the 
number distribution of sources within this period range. 
Assuming the steady formation of sources over time period 
T, the luminosity distribution can be written in the form 

dN t(P) , . dP 
oc AJ-n(P)- 



dL^ 



T 



dL^ 



(Note that the first factor t/T would not appear for a 
LF constructed from a sample of sources produced in an 
instantaneous star formation burst; however, this is the 
case neither for the galactic bulge nor for elliptical galaxies 
we consider here.) 

Clearly, the shape of LF will be determined by the 
dependence of all three factors in this expression on the X- 
ray luminosity. Consider these factors separately for bright 
(with giant donors) and dim (with main-sequence donors) 
LMXBs. 

At luminosities below L x « 10 37 erg/s (M < 10 -9 
M0/year), the majority of donor stars in our LMXB sam- 
ple are main sequence stars with long life times (>Gyr). 
The stable mass transfer episode is also very long, so the 
factor t/T should not strongly depend on L x . The period 
distribution n(P) is generally determined by the evolu- 
tionary history of binaries which includes as a minimum a 
supernova explosion to produce the neutron star, a com- 
mon envelope phase, etc. (or involves dynamical processes 
in dense stellar clusters), and can be found for example 
from population synthesis studies. Let us put it in the 
form dN/d log P and leave it as it is for a while. 

The luminosity of a LMXB due to magnetic stellar 
wind braking depends on the orbital period as a power law 
L x oc P 2 - 5 (see e .g. |Iben fc Tutukov||l984[) or L x oc P 33 
( |Patterson|[T984l ), or L x oc P 3 8 ( |Cote fc Pylyser||l989| ) so 
written in the form d log Pj d log L x this factor is indepen- 
dent of L T . So we find 



n(L x ) oc 



1 dN dlogP 1 dN 



L x d log P d log L x L x d log P 



This would match the observed shape of the LF 
of dim LMXBs if the factor dN/d\ogP ~ const (Z^). 
Unfortunately, it is difficult to construct this distribu- 
tion from observations due to many selection effects. 
Population synthesis calculations (e.g. Fragos et al.| [2008) 
produce a variable distribution, which only roughly can be 
considered as constant. Observations suggest that initial 
orbital periods of binary stars do follow this dependence, 
i.e. dN/dlogP =const (Popova, Tutukov, & Yungelson 



1982), and we shall assume this to hold approximately at 



all stages of binary evolution. 

LMXBs with luminosities higher than L x > 10 37 5 
erg/s have predominantly giants companions (see above). 
The duration r of the mass transfer (and hence the accre- 
tion stage) in these systems is significantly smaller than 
that for systems with main sequence donors, and to a large 
extent is determined by the mass transfer rate (see e.g. 
Webbink, Rappaport, & Savonije 19831). Let us assume 



that the luminosity of LMXB is a powerlaw function of 
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its orbital period L x oc P a (which is true for the magnetic 
stellar wind braking dominated regime, a ~ 2.5 — 3.8, 
see above, and also holds when nuclear evolution of the 
giant donor is responsible for its Roche lobe overflow, 
a ~ 1, 



see 



Webbink, Rappaport, & Savonije 1983). We 



also adopted that the duration of the bright LMXB stage 
is inversely proportional to its luminosity r oc L~ x . This 
is approximately correct because the total mass, which 
might be accreted from a late-type giant donor stars does 



not vary much from system to system (see e.g. Webbink, 
Rappaport, fe Savonije|198"3 ) . Combining these we obtain 



dN 
dL x 



1 dP 1 d\ogP 



oc 



L X P dL x L\ d log L x 



Thus almost independently of the mechanism which 
drives Roche lobe overflow in luminous LMXBs (with late- 
type giants) we obtain a slope of the LF similar to the 
observed one ( |Gilfanov|2004"l |Kim fc Fabbiano|2004| . The 
condition for this to be correct is that the donor stars 
in these LMXBs are short-living giants. We can try to 
estimate the longest lifetimes of such giants assuming that 
the maximum mass, which might be accreted from them 
is ^ O.6M (e.g. |Webbink, Rappaport, fe Savonije||1983 ) 
and their duty cycle is close to unity. In this case r g i ants ^ 
O.6M /(2 x Kr 9 M /year) - 3 x 10 s years. This means 
that the list of bright LMXBs in the Galaxy should alter 
after approximately 30 Myrs. 

We would like to note here that the boundary between 
different types of donors at luminosities above and below 
L x ~ 10 37 erg/sec was also previously noted in numerical 



simulations of LMXB populations in Fragos et al. ( 2008 ) ; 



Kim et al. (2009), but the slopes of the LMXB LF above 



and below this luminosity was found to be not so different. 
The inability of particular numerical simulations to repro- 
duce the observed break in the X-ray luminosity function 
of the LMXB population might indicate that we do not 
(yet) understand all details of the physical processes (e.g. 
the mode of the angular momentum loss, the common en- 
velope phase, etc.) that shape the LMXB LF. We plan to 
perform more detailed calculations of the LMXB popula- 
tion in our future work. 
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